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Abstract. 

The electric dipole moments of closed-shell atoms are sensitive to the parity and time- 
• reversal violating phenomena in the nucleus. The nuclear Schiff moment is one such property, 

" ^ it arises from the parity and time reversal violating quark-quark interactions and the quark- 

chromo electric dipole moments. We calculate the electric dipole moment of atomic 199 Hg 
^ j-H arising from the nuclear Schiff moment using the relativistic coupled-cluster theory. This is 

O |' the most accurate calculation of the quantity to date. Our calculations in combination with 

the experiment data provide important insights to the P and T violating coupling constants 
at the elementary particle level. In addition, a new limit on the tensor-pseudo tensor induced 
atomic EDM, calculated using the relativistic coupled-cluster theory is also presented. 

in 
o 

^ , 1. Introduction 

The origin of parity violation can be explained within the frame work of the Standard Model 
(SM) of particle physics through the weak interactions, but there is no clear understanding 
of the origin of time-reversal violation in nature. The presence of a non-zero electric dipole 
moment of a non-degenerate physical system is a direct signature of parity (P) and time- 
reversal (T) symmetry violations [HO [3]. In this paper we present a study of atomic EDMs of 
closed-shell atoms which arise mainly from the tensor-pseudo tensor (T-PT) electron-nuclear 
interactions and the nuclear Schiff moment (NSM). At the elementary particle level, the origin 
of closed-shell atomic EDMs is attributed to the P and T violating electron-quark interactions 
and quark-quark interactions which are predicted by the lepto-quark models, SUSY, etc [3]. 
The limits on the T-PT coupling constant (Cy) has been obtained from the comparison of the 
most recent experimental result of 199 Hg atomic EDM [5j [6] , |d( 199 Hg)| < 2.1 x 10 _28 ecm and 
the enhancement factor (ratio of atomic EDM to the coupling constant of the interaction in 
question) calculated in [7]. The calculation uses Coupled-perturbed Hartree-Fock (CPHF)[8j[9] 
theory and the calculated value of the atomic EDM is a!( 199 Hg) = —6.0 x lCP 22 C;r<7jvem 



o 



c3 



combining with experimental result, a limit on Ct can be derived, 



C T < 3.26 x 10~ 9 a N 

where on is the nuclear spin. A non-zero value of Ct would imply physics beyond the Standard 
Model. Improving the accuracy of the calculations would lead to a more stringent bound on 
Ct- The CPHF theory accounts for two-particle two-hole type of electron correlations to all 
orders in perturbation. In addition to this, for accurate atomic calculations [lOj. it is important 
to include other correlation effects which are absent in the CPHF theory. Coupled-cluster (CC) 
theory [111 [T2] is one of the many-body methods which is non-perturbative and incorporates all 
forms of correlation effects \13 \ 114 1 [15]. CC theory with all levels of excitations is equivalent to 
all order Many-body Perturbation Theory. It has been successfully applied for calculation of a 
variety of atomic properties |14 } 116 } \T7\. It has also been applied to atomic Hg for calculating 
transition energies [13j . 

The NSM(denoted by operator S) can originate from the nucleon-nucleon interactions or a 
nucleon EDM. At the elementary particle level it can arise from the interaction between the 
quarks and the chromo EDM of the quarks. The coupling constants associated with these 
interactions can be predicted by Multi-higgs, SUSY [H [TH] etc. The dependence of the T-PT 
and NSM interactions on the nuclear spin makes closed-shell atoms, in particular, those having 
non-zero nuclear spin the best candidates to look for EDMs sensitive to the nuclear sector. For 
199 Hg, the EDM induced by the NSM is calculated and parameterized in terms of the Schiff 
moment operator S. The most recent calculation [19, 20J gives 



( ^ 



d Hg = -2.8 x 10- ^ j ecm. 

According to the recent work [2T] the expression for the P and T - odd nuclear potential arising 
from the nuclear Schiff moment 

HR) = -3^p(R), (i) 

where B = J R 4 p(R)dR and R is the electron coordinate and p(R) is the nuclear density. The 
interaction Hamiltonian of electrons with this potential is 

H S u = 3e^p(R) (2) 

where e is the charge of the electron. If p(R) is the normalized density function, which is 1 for 
Rn — 5 < R < Rn + 5, where Rn is the nuclear radius and 5 is a small number. The single 
electron matrix elements of H^m are 



oo 

(ipks 1/2 \HsM\^m Pl/2 ) = 3Se(-l/3) J Pk Sl/2 (r)P mpi/2 (r) +Qks 1/2 (r)Q 
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where Pi and Qi are the large and small components of the relativistic Dirac-wavefunctions 
The matrix elements of Hsm between the states \ip m p 1/2 ) and \tpks 1/2 ) is identical to the 
above expression. 



2. Coupled-cluster theory for closed-shell atoms 

The starting point of obtaining the coupled-cluster equations is the relativistic atomic 
Hamiltonian in the Dirac-Coulomb approximation in atomic units (m e = 1, |e| = 1 and h = 1) 



N N , 

H = J2 \ccti ■ Pi + (A - l)c 2 + Vfr(ri)l + E — ' ( 4 ) 
i i<j Ti i 

where c is velocity of light, a and f3 are the Dirac matrices, l/r^ is the Coulomb potential 
energy between two electrons, and V/v( r «) i s the nuclear potential. In the above Hamiltonian, 
the rest mass energy is subtracted from the total energy eigen values. This is the Hamiltonian of 
an atomic system considering only the electrostatic interactions. The single particle equations 
are obtained by approximating the two-electron term in Eq.Q by a central field potential 
Udf(t), known as the Dirac- Fock potential \22\ 123]. then 

N 

F DC = [ C(X i ■ Pi + (A - !)C 2 + VN(n) + t/ DF ( n )] + V es . (5) 

i 

The residual Coulomb interaction 

N 1 

i<j % 3 i 

The non-central (or) correlation effects are included by treating V es as a perturbation. The 
single electron wavefunctions ip a ) satisfy the Schroedinger equation 



cai ■ Pi + (A - l)c 2 + V N {n) + U DF {n)\ \if) a ) = e a \ifj a ) (6) 

and are expressed in the relativistic framework in terms of the two-component form disucssed 
in succeeding sections. The Eq.([6]) is obtained by variational extremization [22]. The perturbed 
equations are obtained by introducing the P and T violating tensor-pseudo tensor interaction 
Hamiltonian in addition to the residual Coulomb interaction. The exact atomic wavefunction 
in CC formalism |^} = e T |<3?o); where |<&o) is the reference state, T = T\ + T2 + • • • + T/v is the 
cluster operator. In our calculations, we use the singles and doubles approximation T = T1+T2, 
the second quantized form of the operators are 

T X = ^ata a tP|$ } 

a,p 

T 2 = E ^ 44°& a » t Si $ o>. (7) 

a,p,b,q 

The many-body Schroedinger equation of the Dirac-Coulomb Hamiltonian for an atomic 
system, in a state |^) in the CC formalism 



He T \<Z> ) = Ee T \$ ) 



Operating from left side by e 



-T 



e - T #e T |$ ) = E\$ Q ). 



Expressing H in normal ordered form, H = Hn + Epp, where Ep>p = ($o\H\$o) — {&o\Hn\$q) 
is the Dirac-Fock energy. Then, 

e~ T (H N + E DF )e T \^) = E\$ ) (8) 

Projecting Eq.([8j) with singly and doubly excited states an d (& r a t\ respectively, the single 
and double excitation cluster amplitude equations are [21] 

«|{(Se T ) c }|$ ) = (9) 
($;||{(^e T ) c }|$ ) = 0. (10) 
2.1. EDM perturbed coupled- cluster equations 

The -f^EDM operator is the general EDM operator. In the present paper, -£Tedm denotes ffsM, 
which is a single-particle operator. Consider the NSM perturbed Schroedinger equation 

H\W)=E\V) (11) 

where H = H-qc + AHedm arid \^f) = e T |$o) = e T( ° )+AT(1> |$q). Taking upto one order in A in 
the exponent, 

He Tm (l + ATW) |$ > = Ee Ti0) (l + AT«) |d> >. 

Substituting for H in the above equation and comparing zeroth and first order terms of A on 
both sides, 

(ito T(0) )|$o)=£e r( > >, (12) 

and 

(tf DC e T(0) T« + H EDM e TW ) |d> ) = Ee Tm T^ ). (13) 
Multiplying Eq. (fP2"]) by on both sides and substitute H^c = Hn + Ep>p 

TMH N e Ti0) \* ) = AiWT« e r( > ). (14) 
and Eq. (|13p can be written as 

(H N e Tm TW + H EDM e TW ) |* > = (aE^^T®) |$ ) (15) 

since T^ ) and commute and AE co „ = E — Edp is the correlation energy. Operate Eq. (|15p 
on both sides by e _T(0> 

(H N T^ + H EDM ) |$ ) = A^T^I^o) (16) 
Operating by e _:r(0) on EqO 

T^H N \^ ) = AW (1) |$o). (17) 
Subtracting Eq. ([16p from Eq. ()17p 

'h n ,T^} |*o) = -#edm|$o), (18) 



where O 



-t(o) 



Oe T (°) and 6 is any operator. The equation for the -Hedm perturbed singles 
and doubles cluster amplitudes can be derived from the basic equation, Eq. ([T5]) by projecting 
on both sides of the equation with singly and doubly excited determinantal states. 
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l*o> 
l*o> 



which are equivalent to 



Further expanding 
linear matrix equations 




-<^|#edm|$o> 
-(*SI^edm|*o) 

•(^|F EDM |$o), 
■(*3|^edm|*o>. 



(19) 



(20) 



these equations can be cast in the form of a system of 



HuT^ + H 12 T 2 {1) 

U 1^0-) i u T">(1) 



-H20, 



(21) 



where Hn, H12, H21, H22 are the sub blocks of the dressed matrix elements and Hiq and 
H20 are the sub blocks of the dressed Hedm matrix elements. 

2.2. Calculation of Atomic Electric Dipole moments 
The EDM of the atom in the state 1$) is 



D 



atom 



(*l*> 



(22) 



where D is the electric dipole operator. Substituting the above expression for |^) and keeping 
terms only of order A, we obtain the expression for EDM (the symbol of contraction appears 

using where the curly brackets refer to normal ordering and the 

expectation value of normal ordered operator between the vacuum states is zero) 



D 



<$o| 



l$o> 



atom 



(*o|*o) 



(23) 



which can also be written as 



D 



atom 



($o|[5r (1) +r (1)t £>]|$ > 



(*o|*o) 



(24) 



where D = e Tm De T(0) . Using the fact that and D operators are odd and operator 
is even under parity, the bra and the ket vectors in the above expression must have the same 
parity. To simplify the calculations, we expand D as [25, 



D 



1 + 



T (o)t + 



T (o) 



2! 



+ 



71=1 ' ^ ' 



T (0) 



(25) 



In Eq. ()25|) . the one-body nature of the electric dipole operator D restricts the maximum possible 
contractions with to just two. Define 



then 



|$i) = r«|$ ) 

($o\D\$i) + (*i|C|*o 



if + t«) |$ ), 



,(*i|D|*o 



^ atom - (*o|*o> " (*o|*o> ■ 1 j 

The last step follows as the two terms are the complex conjugates of each other and give equal 
contributions. Substituting the expanded form of D 



D 



atom 



De 



T (0) 



n=l 



De 



T (0) 



|$o>/(*o|*o> 



00 1 / ±\ n 



n=l 



/<*o|*o) (27) 



The complexity of the above expression can be mitigated by exploiting the fact that not all 
the terms containing the operators contribute to the infinite summation. In this scheme, 
the zeroth order -D a tom consists of terms without operator, first order -D a tom consists of 

terms having one order of and so on. The unlinked terms of the numerator cancel with 

the denominator and only linked terms contribute in the numerator. At the linear level, 



D 



D + DT W + T (0) D 



= D + DT[ Q) + DT£> + T[ u > ' D + T 2 W ' D 
Following are the terms contributing to the EDM expectation value at the linear level 

Datom = (<t>o\DT^ +T^D\<t> ) 

The two terms in Eql28l are complex conjugates of each other. Hence 

D atom = 2($ |T( 1 ) t D|$o). 
Substituting for D the expression for the atomic EDM 



,(o) 



(o) 



,(o)1 



(28) 



D 



atom 



2($ 



T^D + T^DTP + T^Wf + T 2 (1) W 2 (0) + T^DT} 0) 



|$r 



(29) 



3. Results 

3.1. Basis 

The single particle orbitals for all calculations in the subsequent sections were generated using 
the Gaussian basis [T7] set expansion whose salient features are presented in this section. In the 
central field approximation, the solution of the Dirac equation in terms of the four component 
spinors is given by 

1pnKm{r,6,(j)) = r~ l 



where P n K( r ) an d QnK,{ r ) are the large and the small components of the radial wavefunctions 
expanded in terms of the basis sets [26], 

P nK (r) = Y, C K P 9 L Kp {r) 

p 

Qn.(r) = £C*<4(r) 
p 

where the summation over the indexp runs over the number basis functions N, g^ p (r) and g^ p (r) 
correspond to the large and small components and C% p and C^ p are their expansion coefficients 
for each value of k. The functions g^ p {r) are chosen to be the two parameter Gaussian Type 
Orbitals (GTOs) [26\ and large and small components are related by the condition of kinetic 
balance [271 128] . The calculations presented in this paper correspond to the Even Tempered 
(ET) basis set. For the nuclear density we use the two parameter Fermi distribution 

'"M = 1 + e (r-c)/« 

where r is the radial coordinate, a = 2.3/41n3 is the skin thickness parameter and c = 
\/(5 r rms/3 — 7a 2 7r 2 /3), the half-charge radius [29J . In the expression of c, the r rms is the nuclear 
mean-square radius. The quantity po is determined by normalising pAr(r) over a spherical 

volume. The radial grid has the form, r/% = tq ^ k ^^ h — l\,k = 1, ,n p , n p is the total 

number of grid points. 



Table 1. No. of basis functions used to generate the even tempered Dirac-Fock orbitals and 
the corresponding value of «o an d f3 used. The total number of active orbitals are shown in 
the brackets for Active holes. 





Sl/2 


Pl/2 


P3/2 


^3/2 


d-5/2 


/5/2 


fl/2 


97/2 


99/2 


Number of basis 


31 


32 


32 


20 


20 


20 


20 


10 


10 


ao(xlO~ 5 ) 


725 


715 


715 


700 


700 


695 


695 


655 


655 


P 


2.725 


2.715 


2.715 


2.700 


2.700 


2.695 


2.695 


2.655 


2.655 


Active holes (36) 


2 


2 


2 


2 


2 


1 


1 


1 


1 


Active holes (39) 


3 


3 


3 


2 


2 


1 


1 


1 


1 


Active holes (43) 


3 


3 


3 


3 


3 


2 


2 


1 


1 


Active holes (45) 


3 


3 


3 


3 


3 


3 


3 


1 


1 


Active holes (51) 


5 


5 


5 


3 


3 


3 


3 


1 


1 


Active holes (57) 


7 


7 


7 


3 


3 


3 


3 


1 


1 


Active particles 


6 


4 


4 


3 


3 


1 


1 









For the present calculation, the details of the basis set used in generating the orbitals is 
shown in Table. [H Using the above basis, the single particle orbitals have been generated. Next, 
an existing closed-shell coupled-cluster code developed inhouse has been used to generate the 
unperturbed cluster amplitudes [30]. The code for calculating the perturbed cluster amplitudes, 
and subsequently the atomic EDM was developed by our group [25j 124"] . 



Table 2. Individual contributions 



Contributions in units 
of e ao Gf Ct 





-47.830 


T (i^ DT (0) 


0.082 


t^dtP 


14.320 


T (^ DT (0) 


-0.386 


T (^ DT (o) 


-0.059 



Total -33.874 



3.2. Results for Hg EDM induced by the P and T violating T-PT interaction 
Contributions from each of the terms in Eq. (|29p are shown in Tablej2] for the basis with 57 
active orbitals. The final result can be expressed in units of e-m 

£> Hg = -1.125 x W- 22 C T ema N . 

The Dirac-Fock contribution is —2.45 x 10~ 22 emCT<r/v- It can be noticed from the Table[2]that 

the largest contribution is from D. This can be explained from the fact that it includes the 
Dirac-Fock effect, which is the most dominant contribution. The trend shown by the individual 
contributions in Tabled is related to the fact that the cluster amplitudes are smaller in 
magnitude compared to the cluster amplitudes. In addition, the T± amplitudes are 
larger in magnitude compared to the T% amplitudes again due to the presence of the Dirac- 
Fock contribution in For example, from the above table, we see that the contribution of 
DT^ is larger than that of DT^ . Similarly, the contribution of DT^ 1 is greater 

than that of DT^ . These arguments are equally valid for the atomic EDM induced by 
the nuclear Schiff moment shown in the next section, which follows the same trend as above. It 

is also interesting to note that in both cases contribution of the term T^' D is approximately 

3 times larger in magnitude than that of the term DT^\ which is the second largest 
contribution. 

3. 3. Results for Hg EDM induced by the P and T violating Nuclear Schiff moment 

The 199 Hg atomic EDM induced by the nuclear Schiff moment is calculated with the same 
basis given in the previous section. The methods of generating the perturbed and the 
unperturbed cluster amplitudes are the same as described earlier. The Dirac-Fock contribution 
is £>H g = —0.546 x 10 5 e a,Q S = —0.390 x 10 _17 ecm5(e/m 3 ) -1 . Contributions from each of 

the terms in Eqj29] is shown in TablelHJ Again, the leading contribution is from if 1 D and 

it can also be seen that the contribution of D is approximately 3 times in magnitude of 

the contribution of DT^, which is also true for the results of the T-PT induced EDM. 



Table 3. Individual contributions 



Contributions in atomic 
units of S e xqq 



T[ l)] D -0.177 xlO 5 

T Q) DT (0) omo xlQ 3 

T[ l)] DT^ ] 0.525 xlO 4 

T 2 (1) Z>T 2 (0) -14.258 xlO 1 

T 2 (l)1 DTP 0.252 xlO" 4 

Total -0.126 xlO 5 



Our result is not in agreement with Dzuba et.al [19]. They have used CI + MBPT method 
for the generation of the orbitals. We have compared the Schiff moment interaction and the 
electric dipole matrix elements of the 6si/ 2 an d the core pi/ 2 orbitals for 199 Hg with the results 
obtained by the authors of Dzuba et al. at the Dirac-Fock level and found that the agreement 
was very good. This suggests that the discrepancy in our results could be majorly due to the 
different choices of the virtual orbitals in the two calculations. 
The final result is in the units where S is expressed in e/m 3 

= 2 x 10~ 23 S 
D Ug = -0.126 x 10° x 2 x ——recm 

0.529 z ejm 6 

S 

D H = -0.0901 x 10" 17 ecm 



e/m 3 



3.4- Summary of the EDM results of 99 Hg 

The interaction Hamiltonian for the T-PT and the NSM are both dependent on the nuclear 
density /9/v(r) and hence their matrix elements are sensitive to the S\/2 an d the P\i2 orbitals, 
which have a non-zero probability density inside the nuclear radius. In addition, the NSM 
interaction Hamiltonian is proportional to the electron coordinate R. In addition, the matrix 
elements of the atomic EDM contain the electric dipole operator D. Hence the atomic EDMs 
arising from these interactions require the single particle orbitals to be very accurate at all 
radial ranges. The trend followed by the T-PT and NSM induced EDM shows that the atomic 
EDM is very sensitive to the inclusion of Si/ 2 and the p\/2 virtuals. On the other hand, there 
is not much variation in the polarizability as it is more sensitive to orbitals of higher orbital 
angular momenta. We have also performed linear CCEDM calculations without including any 
cluster amplitudes. This in other words, amounts to the linear CCEDM calculation with 
the cluster amplitudes generated with the approximations Hn ~ Hn and Hedm ~ -Hedm hi 
Eq. (|20|) , Hence all the correlation effects from amplitudes are omitted. The results of this 
calculation for a basis of 39 active orbitals are discussed. The Dirac-Fock contribution for this 
basis is -DHg = —3.17 x 10~ 22 em. The total contribution is -DHg = —5.77 x 10 _22 em, whereas 
the CPHF calculation gives Z?Hg = —4.64 x 10 _22 em. Similar comparison can be made for a 
basis of 57 active orbitals, for which the Dirac-Fock contribution is -DHg = —2.45 x 10 _22 em, 



Table 4. Summary of our results. 



Basis size 


In units of 
ru ujeni 


In units of 

i r\ — 17 , , , ,-, . , 5" 
iu e cm —r — t 


36 


-2.70 


1.83 


39 


-2.97 


-0.23 


42 


-2.21 


-0.15 


45 


-1.41 


-0.11 


51 


-1.40 


-0.11 


57 


-1.13 


-0.09 



Table 5. P and T violating parameters 



P and T violating parameter at the Limits obtained from 
elementary particle level linear CCEDM 

Vnp 
9irNN 
#QCD 

e(d u - d d ) 



the bare-Coulomb calculation without amplitudes gives Z?Hg = —6.60 x 10 _22 em, the 
CPHF calculation gives Dn g = —5.61 x 10 _22 em. This comparison helps in understanding 
the interplay between the various many-body effects : the the bare-Coulomb which contains all 
kind of correlation effects to one order in Coulomb interaction;the CPHF theory which contains 
the Coulomb effects only of the two particle-two hole kind; the linear CCEDM which subsumes 
the effects of the bare Coulomb, the CPHF and more. 

Limits on the P and T violating parameters at the elementary particle level can be derived 
from the tensor-pseudo tensor coupling constant (Ct) and the nuclear Schiff moment (S) [2"5"||2"4"]. 
The results of various P and T violating parameters are shown in TableEl A detailed discussion 
on the derivation of these parameters has been presented in |25j . Using these important 
parameters, the observable EDMs can be connected to the P and T violating coupling constants 
predicted by the underlying elementary particle theories. The bounds on these coupling 
constants can further be used to constrain various models of elementary particle physics. 
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